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ABSTRACT 

I investigate non-perturbative aspects of zero-dimensional matrix models. Sub- 
tleties in the large- limit of the semiclassical picture are pointed out. The tun- 
neling of eigenvalues is seen to correspond to a chaotic sequence of recursion coef- 
ficients determining the orthogonal polynomials. 
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Introduction. 

Over the last four years, we have learned how to model two-dimensional euclidean 
quantum gravity (with topological fluctuations) by hermitean matrix models [1], 
in the so-called double-scahng limit [2]. The latter entails sending the size N of 
the matrix M to infinity while the couplings g in the matrix potential V approach 
a critical value gc in such a way that some combination of N and (g—gc) is kept 
constant. In the pure gravity case, one has 

V{M) = iM^ + gM^ , Qc^-^ . (1) 

Interestingly, the critical matrix potential retains a local minimum but is un- 
bounded from below, and analytic continuation is to be performed into a region 
where the partition function is ill-defined. The result is either singular or develops 
an imaginary part, associated with the instability towards the tunneling escape of 
individual matrix eigenvalues [3]. 

The tunneling phenomenon is most transparent in a semiclassical treatment of 
large- matrix models. As a first step, the saddle-point analysis was pioneered in 
ref. [4]. To work with a well-defined theory, I will bound the potential (1) from 

below by changing its 'large-M' behavior through a modification 6V{M) = eM^, 
with a small, positive regulator e. This turns out to change substantially the large- 
N phase structure of the model [5-8]; in particular, the critical line gc = gs{^) 
becomes metastable and invisible at finite A^. Nevertheless, tunneling may now be 
studied directly, without the need for analytic continuation. 

Instrumental to the success of the random matrix model has been the tech- 
nology of orthogonal polynomials [9,10]. The saddle-point results, however, have 
not yet been completely understood in this framework. Only in the case of degen- 
erate potential wells, a precise relation between orthogonal polynomial recursion 
coefficients Rj^ on one side and the (multi-band) classical eigenvalue density on the 
other side has been established [11,7,8]. 



1 



The more general connection between the two approaches involves an inter- 
pretation of eigenvalue tunneling in terms of orthogonal polynomials. This will be 
addressed in the third part of my talk, where I shall present a surprising resolution 
of a puzzle mentioned earlier [7,8,12]. The second part, following this introduction, 
will outline the semiclassical approach. Along the way, I will relate N—oo saddle 
points to airfoils, and interpret eigenvalue tunneling in terms of the thawing of 
a frozen Dyson gas. A short summary shall close my talk, which reviews some 
research of mine conducted over the past three years and published for the most 
part in refs. [8,12,13]. 

The Semiclassical Approach. 

I begin by formulating a collective field theory for the zero-dimensional hermitean 
one-matrix model at finite N. My starting point is the partition function 

Zn oc y"Ae-^*^^(^) (2) 

for an NxN random hermitean matrix ensemble, in a potential V. Upon diago- 
nalization M = diag{xi) this reduces to 



Ylfdxi exp|-iV^y(a;i) +^ln(xi-Xj)4 , 



(3) 



the partition function of a two-dimensional Coulomb gas of charges restricted to 
a line, in an external potential V{x) at temperature /3^^=l/iV^. I like to change 
variables from the matrix eigenvalues Xi to their density distribution 

N 

P(^) - ^ I] - ■ (4) 

i=l 

More precisely, I insert 



1 = [Vpll6{p{x)-j^j:^5{x-xi)) 

= JJ VpVX expji J dxX{x)[p{x) - j^J2iS{x-Xi)]^ 



(5) 
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into (3) and express the action in terms of the density, 
S%[p\ ^ N'' jdxp{x)V{x)-\N^ jjdxdyp{x)f{x-y)p{y) + \Nf{Q) . (6) 

The self-interaction had to be regulated by replacing Inz^ — > f{z) in eq. (3), 
choosing some suitable, i.e. symmetric and bounded, function /. 

Following ref. [14] I am able to perform the integration over Xi, 



Zn = ff VpVX e-^^[^]+*/^^ Yl Idx 



(7) 



VpVX 6"^^!'']+^/^^ 



dx e N 



and arrive at an exact effective action 



Sn[p,X\ = S%\p\ - i J dxX{x)p{x) -N\n j dx 



(8) 



which is not only nonlocal in the two real fields p and A but also non-polynomial 
in the latter. Interestingly, the constant mode of A can be integrated out exactly 
to yield the constraint S{Jp—l) that was apparent already from the definition (4). 
However, I shall keep those modes in the measure for the time being. In principle, 
another constraint arises from the positivity of p. Perturbation theory about a 
strictly positive p, however, is insensitive to this restriction, and I will therefore 
ignore it in the following. 

My goal is to initiate a systematic semiclassical analysis of this peculiar one- 
dimensional field theory. To leading order in h we must determine the saddle-point 
configurations (p. A), where the action (8) is stationary. The first variations yield 

= W(x)-7V2 jdyf{x-y)p{y)-iX{x) 



(9) 

= p{x) - e-y^^''^ / Jdye-y^^y^ 
where p comes out to be properly normalized, J p — 1. The second equation 



3 



determines A up to a constant, 

A(x) = iN In p{x) + Xo , (10) 

a result which may be inserted into the first equation* Differentiating with respect 
to X to remove constants and deregulating the Coulomb repulsion I get^ 

where ^ denotes Cauchy's principal value of the integral. The r.h.s. is known as the 
Hilbert transform H (of p) which has been thoroughly investigated [16]. Together 
with normalization and positivity, this equation describes the classical eigenvalue 
density for any finite N. It is noteworthy that (11) is not homogeneous in so its 
solution cannot be, either. At N — oo the equation (11) has been widely studied and 
solved [17,18,4], and it was learned [5,12] that a unique solution extends to < oo. 
Unfortunately, the equation is not easily solved for finite A^. Even in the large- A^ 
limit some care is required, as shown by the following. For instance, should I drop 
the p'/p-tcrm since it is down by -^? A little inspection reveals that such a step 
is in general not consistent with the asymptotic large- behavior of the equation, 
which demands a - fall-off for the l.h.s. due to the normalization of p. In fact, the 
formal 'solution' 

Pe{x) - e^p{-N[V{x) - JH^-yf Peiv)]} '"fe" x^^e-^^^^) (12) 

of equation (11) shows that the ^ term dominates the l.h.s. of eq. (11) for |x| 1 
(unless V ~ Inx^), so that p^ = 0{e^^) asymptotically. Note that I have attached 
a subscript e = to the solution to indicate that it depends on the value of A". 
One can read off that the effective potential seen by eigenvalues |x|^l approaches 

★ The first equation subsequently fixes the constant Aq for a given solution p. 
t See also ref. [15]. 
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V{x) — 21n|x|. The situation is different, however, near the minimum of the 
potential where V' dominates the l.h.s. of eq. (11) and most of p is concentrated? 
When p is (9(1) the 4-term may safely be neglected or treated as a perturbation 
in (11), leading to 

by simply inverting the Hilbert transform in equation (11). The result is a mod- 
ulation of Wigner's semicircle distribution. When — > oo, the crossover regions 
between the interior and exterior of the 'Fermi sea' shrink to points x—a, b, and 
the saddle-point equation (11) turns into [4] 





W'{x) = /-^Po(l/) = T^Xcipo] fora;e[a,6] 



a 



(14) 

po{x) — for a; ^ [a, 6] 

This relation is known as the airfoil equation and determines the vorticity po related 
to a given velocity field V' along the airfoil {b—a is the span of the wings) [17]. 
In our case the location of a and b is determined from the normalization of po- 
Equation (14) is solved by inverting the finite Hilbert transform T, 



^/{b — x){x — a) 



TT tt"^ J X — y 



■ (15) 



a 



Again, I have assumed a single-well potential V, so that the support of po is a 
single, connected interval [a,b]. 

It is worthwhile to give the form of the saddle-point action. Employing eq. (10) 



I The case of multiple local minima will be discussed in a while. 
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— — X 

as well as In J e"-^^ — —j^Xq I find, in agreement with ref. [15], that {p—p^) 

Sn[p,X] = N'^ j p{x)V{x)-N^ jj p{x)\n\x^\p{y) + N j p{x)\np{x) 

= ^ j p{x){V{x)-\u\x\) + ^V{Q) + ^ Jp{x)\np(x) + ^\np(0) 

(16) 

where I dropped a (singular) term ^Nf(0) and made use of the saddle-point equa- 
tion (11). The last integral permits an interpretation as the entropy of the distri- 
bution p. Taking the naive large- limit, I obtain 

SoolPo] = ¥ Jdxpo{x)V{x)-N^j , (17) 



with 

7 = J dy hi\x — y\ po{y) — ^V{x) — constant for x e [a, b] (18) 

being the chemical potential (or Lagrange multiplier enforcing J po — 1). 

At this point, I would like to drive home an essential point of my talk. The 
non-trivial A'^-dependence of the classical backgound pe implies that the classical 
limit, p Pe, differs from the low-temperature limit, — > oo, because of cor- 
rections coming from the integration measure. As a consequence, the semiclassical 
loop expansion will not be identical to the topological expansion; rather, a dou- 
ble expansion arises. Beyond this, non-perturbative (in -^) contributions appear 
already at tree level (in K). It is therefore by no means clear that interchanging 
limits, by first taking — > oo in equation (11) and then solving equation (14) 
to obtain po, provides the large- A?^ limit, lim£_»o Pe, of a proper solution to equa- 
tion (11), as my notation suggests. In fact, I will now demonstrate that in general 
a solution po does not correspond to a finite- A?^ saddle point p^. 

Non-perturbative effects become tangible when two or more potential wells 
compete for eigenvalues. For even potentials, the full complexity of the problem 
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appears first in the triple-well potential 



V{x) = Ix^ + gx'^ + ex^ , e>0 , (19) 

since the x <-> —x symmetry trivializes the double-well case. I like to fix e to some 
small value and probe the phase diagram by decreasing the quartic coupling g 
along the negative axis. For g < = — ^3e/2 the potential develops three well- 
separated minima which become degenerate &t g = g@ = —\ple. For sufficiently 
negative g the 'Fermi sea' must, therefore, consist of three or two disconnected 
oceans, called arcs or bands. At A?" ^ oo, the eigenvalue density should then be 
supported on one, two, or three disjoint intervals. As shown in refs. [6,7,8], the 
general solution of equation (14) reads 



Po 



(x) = f {n^ - x^) ^J{4a^ - x^) (4^2 - x^) (4c2 - x^) , (20) 



describing a positive three-band density with < 2a < n < 2b < 2c and sup- 
port on [—2c, —2b] U [—2a, 2a] U [26, 2c]. Normafization imposes three conditions on 
{n, a, b, c} which leaves a one-parameter family of solutions. A convenient parame- 
ter to label these solutions is the difference A7 = ^o—lc of the chemical potentials 7 
for the outer and the central bands. The chemical potential is nothing but the in- 
tegration constant appearing when integrating equation (14). Since it may take 
different constant values for two eigenvalues Xc G [—2a, 2a] and Xq G [26, 2c] from 
two different bands, the difference 

A7 = (dy In 1^ pod/) - \ [V{xo)-V{x,)\ (21) 

J \Xc-^\ 

is a genuine property of the solution po. In this situation one must replace 

bands ^^ij ^^^^ bands 

in equation (17). Extremal values of A7 occur when the number of bands decreases, 
for 6 — > c or a — > 0, and the density (20) becomes unique. However, those solutions 
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do not exist everywhere in the (e, g) plane. For sufficiently small e, the range of one-, 
two- and three-band solutions is given by the sequence < g@ < g2 < 9* < 91 < 0. 
Here, one-band densities arise for g > gs, ^w;o-band distributions occur for g < 92, 
and three-hand solutions appear for g < gi. Hence, equation (14) admits a unique 
solution only for g > gi. The overlapping regions above indicate a coexistence 
of multiple-band densities elsewhere. Even more astonishing is the discovery that 
the members pQ{x,A'y) of such a family are not degenerate in free energy [8]. 
This contradiction in terms is resolved by noticing that infinitesimal variations 
within the family, which correspond to the tunneling of individual eigenvalues, 
are actually singular at the band edges. One might say that it requires a finite 
variation to move an individual eigenvalue to another well although it is only 
an 0{jf) effect. Since -^^po is not square-integrable the inclusion of this mode 
among the density fluctuations in debatable. As the chemical potential drives the 
tunneling of eigenvalues, one should expect the minimal action to belong to the 
unique family member with A7 = 0, which I call dominant. This is indeed what 
happens and can be checked numerically [8]. 

Of course, for < 00 there is always a unique saddle-point density pe, with a 
unique limit as e ^ 0, because the strict positivity of the distribution implies that 
A7 = all along. Hence, we have 



With hindsight it is clear that the 'sub-dominant' members of a three-band family 
could only appear because I took the limit A" — > 00 prematurely by going from 
equation (11) to (14). 'Physically' speaking, the freezing of the Dyson gas of 
eigenvalues at zero temperature entirely suppresses any tunneling and permits 
those ficticious saddle-point distributions. At finite temperature, the Dyson solid 
melts at the edges, and tunneling, although exponentially small, destabilizes all but 
the dominant solutions. The 'entropy term' in equation (11) plays the crucial role. 

* 9 = 93{^) is the critical BPIZ line where the double scaling limit is to be taken [4]. 




(23) 
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However, I may still employ the incomplete large- saddle-point equation (14), if 
it is complemented by the additional global requirement A7 = 0. 

In this hght the correct phase diagram for the potential (19) looks as follows. 
Again taking e small but fixed, I find a new sequence g" < g@ < g' < 0, where a 
small region of three-band dominance around the degenerate point {g = g@) sepa- 
rates two-band from one-band densities. Inspecting the shape of the potential at 
the transition values g' and g" , one learns that essentially the number of degen- 
erate absolute minima determines the number of eigenvalue bands. Interestingly, 
the jump from one to three well-separated bands is smoothed out by an inter- 
polating three-band region occuring when non-degenerate wells are in some sense 
comparable and can compete for eigenvalues. 



Orthogonal Polynomials. 

The standard approach to matrix model calculations, and so far the only one 
capable of producing the topological expansion and the double-scaling limit, is the 
method of orthogonal polynomials [9] . Its starting observation is that the van der 
Monde determinant in 



N 



■1=1 



n 

Kj 



(24) 



can be rewritten as 





1 


1 


1 






Pq{x2) . 


Pq{xn) 




Xl 


X2 


XN 




Pi{xi) 


Pl{x2) . 


Pi{xn) 




xl 


™2 






P2{xi) 


P2{X2) . 


■ P2{xn) 




^1 


^2 


■ -^N 




Pn-i{xi) 


Pn -i{x2) ■ 


. Pn-i{xn) 
(25) 



n 

Kj 



with monic polynomials Pk{x) — x^ + (lower order). If one cleverly chooses the P^ 
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to be mutually orthogonal with respect to the measure e -^^(^), 

hkSki = Jdxe-'^^^'^^ Pk{x)Pi{x) , (26) 

the change of basis — > Pk{x) in equation (24) exactly diagonalizes the Coulomb 
interaction and trivializes the free energy to 

Fn = -^]n[NlhQhih2...hN-i] . (27) 

The construction of the polynomials simplifies thanks to the classic recursion rela- 
tion 

Pk+i(x) ^ xPk(x)-RkPk-i(x) , = >0 , (28) 

n-k-i 

so that is it sufficient to compute the norms hk, starting with the initial condition 
Ro = 0. 

The Stieltjes method of constructing the orthogonal polynomials consists of 
iterating eq. (28) and the norm computation, (26): 

Pk{x; Ro, . . . , Rk-i) — > hk — > Rk — > Pk+i{x; Rq, . . . , Rk) ■ (29) 

Unfortunately, it is quite inappropriate for numerical analysis. However, the special 
form e~^^ of the measure allows for a finite recursion relation among the Rk 
themselves, the 'string equation' [9] 

jf = Rk\^ + 4:g{Rk-i + Rk + Rk+i) + ^^{Rk-i + Rk + Rk+i)'^ 

/ (30) 

+ 6e(i?/c_2-Rfc-l — Rk-lRk+l + Rk+lRk+2) I , 

displayed here for the potential (19). After solving for Rk — Rk{Rk-i, ■ ■ ■ , Rk-A) 
one still needs the initial values 

besides R-i = = Rq to start the iteration. It turns out that this procedure 
is numerically tractable, but Rk becomes increasingly sensitive to the initial con- 
ditions for growing k or N [12]. Trivial but instructive is the exactly solvable 
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example of purely quadratic potential, i.e. g = e = 0, or V = ^x'^. In this case one 
simply rediscovers the Hermite polynomials from Ri^ = k/N, as appropriate for 
the harmonic oscillator. 

Let us now investigate the large- A?" or planar limit, in order to connect up 
with the semiclassical results. Here, I have to rely on an assumption, namely a 
continuum approach of the Rj^ needs to be postulated. Writing 

^ = e e [0, 1] , Rk = niO , e = i , (32) 

the simplest ansatz 

Rk+i-Rk = 0{j^) as N^oo (33) 

implies re(^) a smooth positive function with r(0) = 0. Numerical studies 

show that this behavior indeed occurs whenever a single potential well is clearly 
dominant [12]. This coincides with the one-band regime of our potential.* The 
string equation is then dramatically simplified to the algebraic relation [9] 

^ = r{l + 12gr + 60er^} . (34) 

The condition r(0) = selects a unique branch of r which monotonically reaches 
^ = 1 provided g > gsie), the BPIZ critical line. Finally, the free energy Fqo is 
obtained by naively taking the N ^ oo limit of eq. (27), 



N 



Fn = -^ln7V!-^ln/io-^ J](l-#)lnit:fc 

k=l 

1 



to be compared to -^Soo from equation (17). 



(35) 



★ excluding regions of ill-separated minima. 
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However, the continuity assumption (33) is clearly violated when g < g@, 
because an estimate of Ri from equation (31) reveals that it must jump from 
Oljf) to the square of the location of the outer potential minima, ~ — ^ = 
when g drops below the degenerate point, g — g@- Such a behavior is known from 
studies of double- well potentials [11,7], where an alternating sequence 

,1 , , X I t^^HE) for k even , , 

^O(^) . rM)^[J^ ^^^^^^^ (36) 

is observed. Under this modified assumption, the string equation (30) turns into 
two coupled cubic equations for r^°) and Their graphical solution exhibits 

r^^\C,) as an increasing function starting from r(°^(0) =0, and r^^\$,) as decreasing 
from A^\0) 



^2 
min 



Like for the saddle-point method, it is not evident which large- A^" assumption 
is to be chosen for given values of e and g. It can be shown, however, that the 
series of \^Rk approximates the sequence of consecutive eigenvalues Xk eventually 
building up to the distribution Pe{x) in the limit e ^ 0. Hence, one- and two-band 
regions in the phase diagram must correspond to the continuum behavior of eqs. 
(33) and (36), respectively? An unsettling gap remains, however, in the three-band 
dominated buffer zone. How can the one-branch ansatz (33) merge with the two- 
branch ansatz (36) when A^\0) — A^\0) ~ x'^in — "^"^ '^^^ resolution of this 
paradox was discovered through a numerical analysis, which uncovered a ^ interval 
with seemingly chaotic recursion coefficients Rk, interpolating between one 
or two branches for ^ < ^ and a single branch for ^ > ^(> 1). The three-band 
region g" < g < g' coincides with ^ < 1, i.e. the onset of the unpredictable 
behavior creeping into the relevant ^ interval [0, 1] . Only exactly at the degenerate 
point, g = g@, is a simple three-branch solution realized [8]. This picture has been 
confirmed by several groups [19,20]. 



t Actually, the number of branches r^'^ reaching the point ^ = 1 is relevant. 
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Conclusions. 

I have reviewed some non-perturbative aspects of hermitean random matrix mod- 
els, with an emphasis on the distribution of eigenvalues among several potential 
wells. It turned out that the classical, quasi-continuous, large- A?" eigenvalue distri- 
bution depends non-perturbatively on so that the semiclassical loop expansion 
must be distinguished from the standard topological (or string loop) expansion. 

For multiple- well matrix potentials, interchanging the limits N^oo and (S'jv— *-extremum| 
is dangerous due to the N<oo equilibration between different Dyson gas compo- 
nents, mediated by eigenvalue tunneling. As a consequence, the point of critical 
coupling becomes metastable for pure gravity when V gets bounded from below. 

From numerical simulations of the sequence of recursion coefficients Rk for 
the orthogonal polynomials, 1 conjecture that their behavior is characterized by 
the critical points of the potential function V itself. The consecutive equilibrium 
deposition of eigenvalues on the real line is suggested to be sensitive to features 
of V at increasing values. For degenerate absolute minima one observes a quasi- 
periodic series of Rk, whereas non- degenerate minima produce chaotic behavior! 
The unpredictability of the coefficients reflects the competition of incommensurate 
potential wells for eigenvalues. 

I hope to have demonstrated that the continuum limit of matrix models is more 
complicated than imagined originally. In view of this it would be very beneficial 
to understand their critical properties in the semiclassical description. An attempt 
to push the latter beyond the classical hmit is currently in progress (see [13]). 

Acknowledgements: 1 am grateful to the organizers for the charming and stimu- 
lating atmosphere of the Symposium, and especially K. Behrndt for his personal 
efforts to make it all work. 
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